Proper network randomization is key to assessing social balance.

Social ties, either positive or negative, lead to signed network patterns, the subject of balance theory. For example, strong balance introduces cycles with even numbers of negative edges. The statistical significance of such patterns is routinely assessed by comparisons to null models. Yet, results in signed networks remain controversial. Here, we show that even if a network exhibits strong balance by construction, current null models can fail to identify it. Our results indicate that matching the signed degree preferences of the nodes is a critical step and so is the preservation of network topology in the null model. As a solution, we propose the STP null model, which integrates both constraints within a maximum entropy framework. STP randomization leads to qualitatively different results, with most social networks consistently demonstrating strong balance in three- and four-node patterns. On the basis our results, we present a potential wiring mechanism behind the observed signed patterns and outline further applications of STP randomization.

scDrugPrio: a framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

BACKGROUND: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. METHODS: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. RESULTS: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. CONCLUSIONS: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package ( https://github.com/SDTC-CPMed/scDrugPrio ).

scDrugPrio: A framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

Background: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. Methods: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. Results: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. Conclusion: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package (https://github.com/SDTC-CPMed/scDrugPrio).

Neuroscience Needs Network Science.

The brain is a complex system comprising a myriad of interacting neurons, posing significant challenges in understanding its structure, function, and dynamics. Network science has emerged as a powerful tool for studying such interconnected systems, offering a framework for integrating multiscale data and complexity. To date, network methods have significantly advanced functional imaging studies of the human brain and have facilitated the development of control theory-based applications for directing brain activity. Here, we discuss emerging frontiers for network neuroscience in the brain atlas era, addressing the challenges and opportunities in integrating multiple data streams for understanding the neural transitions from development to healthy function to disease. We underscore the importance of fostering interdisciplinary opportunities through workshops, conferences, and funding initiatives, such as supporting students and postdoctoral fellows with interests in both disciplines. By bringing together the network science and neuroscience communities, we can develop novel network-based methods tailored to neural circuits, paving the way toward a deeper understanding of the brain and its functions, as well as offering new challenges for network science.

Neuroscience needs Network Science.

The brain is a complex system comprising a myriad of interacting elements, posing significant challenges in understanding its structure, function, and dynamics. Network science has emerged as a powerful tool for studying such intricate systems, offering a framework for integrating multiscale data and complexity. Here, we discuss the application of network science in the study of the brain, addressing topics such as network models and metrics, the connectome, and the role of dynamics in neural networks. We explore the challenges and opportunities in integrating multiple data streams for understanding the neural transitions from development to healthy function to disease, and discuss the potential for collaboration between network science and neuroscience communities. We underscore the importance of fostering interdisciplinary opportunities through funding initiatives, workshops, and conferences, as well as supporting students and postdoctoral fellows with interests in both disciplines. By uniting the network science and neuroscience communities, we can develop novel network-based methods tailored to neural circuits, paving the way towards a deeper understanding of the brain and its functions.

A positive statistical benchmark to assess network agreement.

Current computational methods for validating experimental network datasets compare overlap, i.e., shared links, with a reference network using a negative benchmark. However, this fails to quantify the level of agreement between the two networks. To address this, we propose a positive statistical benchmark to determine the maximum possible overlap between networks. Our approach can efficiently generate this benchmark in a maximum entropy framework and provides a way to assess whether the observed overlap is significantly different from the best-case scenario. We introduce a normalized overlap score, Normlap, to enhance comparisons between experimental networks. As an application, we compare molecular and functional networks, resulting in an agreement network of human as well as yeast network datasets. The Normlap score can improve the comparison between experimental networks by providing a computational alternative to network thresholding and validation.

Next-generation large-scale binary protein interaction network for Drosophila melanogaster.

Generating reference maps of interactome networks illuminates genetic studies by providing a protein-centric approach to finding new components of existing pathways, complexes, and processes. We apply state-of-the-art methods to identify binary protein-protein interactions (PPIs) for Drosophila melanogaster. Four all-by-all yeast two-hybrid (Y2H) screens of > 10,000 Drosophila proteins result in the 'FlyBi' dataset of 8723 PPIs among 2939 proteins. Testing subsets of data from FlyBi and previous PPI studies using an orthogonal assay allows for normalization of data quality; subsequent integration of FlyBi and previous data results in an expanded binary Drosophila reference interaction network, DroRI, comprising 17,232 interactions among 6511 proteins. We use FlyBi data to generate an autophagy network, then validate in vivo using autophagy-related assays. The deformed wings (dwg) gene encodes a protein that is both a regulator and a target of autophagy. Altogether, these resources provide a foundation for building new hypotheses regarding protein networks and function.

Assessment of community efforts to advance network-based prediction of protein-protein interactions.

Comprehensive understanding of the human protein-protein interaction (PPI) network, aka the human interactome, can provide important insights into the molecular mechanisms of complex biological processes and diseases. Despite the remarkable experimental efforts undertaken to date to determine the structure of the human interactome, many PPIs remain unmapped. Computational approaches, especially network-based methods, can facilitate the identification of previously uncharacterized PPIs. Many such methods have been proposed. Yet, a systematic evaluation of existing network-based methods in predicting PPIs is still lacking. Here, we report community efforts initiated by the International Network Medicine Consortium to benchmark the ability of 26 representative network-based methods to predict PPIs across six different interactomes of four different organisms: A. thaliana, C. elegans, S. cerevisiae, and H. sapiens. Through extensive computational and experimental validations, we found that advanced similarity-based methods, which leverage the underlying network characteristics of PPIs, show superior performance over other general link prediction methods in the interactomes we considered.

Computational Inference of Synaptic Polarities in Neuronal Networks.

Synaptic polarity, that is, whether synapses are inhibitory (-) or excitatory (+), is challenging to map, despite being a key to understand brain function. Here, synaptic polarity is inferred computationally considering three experimental scenarios, depending on the nature of available input data, using the Caenorhabditis elegans connectome as an example. First, the inputs consist of detailed neurotransmitter (NT) and receptor (R) gene expression, integrated through the connectome model (CM). The CM formulates the problem through a wiring rule network that summarizes how NT-R pairs govern synaptic polarity, and resolves 356 synaptic polarities in addition to the 1752 known polarities. Second, known synaptic polarities are considered as an input, in addition to the NT and R gene expression data, but without wiring rules. These data train the spatial connectome model, which infers the polarity of 81% of the CM-resolved connections at >95$>95$ % precision, while also inferring 147 of the remaining unknown polarities. Last, without known expression or wiring rules, polarities are inferred through a network sign prediction problem. As an illustration of high performance in this case, the generalized CM is introduced. These results address imminent challenges in unveiling large-scale synaptic polarities, an essential step toward more realistic brain models.

Geometry of rare regions behind Griffiths singularities in random quantum magnets.

In many-body systems with quenched disorder, dynamical observables can be singular not only at the critical point, but in an extended region of the paramagnetic phase as well. These Griffiths singularities are due to rare regions, which are locally in the ordered phase and contribute to a large susceptibility. Here, we study the geometrical properties of rare regions in the transverse Ising model with dilution or with random couplings and transverse fields. In diluted models, the rare regions are percolation clusters, while in random models the ground state consists of a set of spin clusters, which are calculated by the strong disorder renormalization method. We consider the so called energy cluster, which has the smallest excitation energy and calculate its mass and linear extension in one-, two- and three-dimensions. Both average quantities are found to grow logarithmically with the linear size of the sample. Consequently, the energy clusters are not compact: for the diluted model they are isotropic and tree-like, while for the random model they are quasi-one-dimensional.

Proteomic Analysis Uncovers Measles Virus Protein C Interaction With p65-iASPP Protein Complex.

Viruses manipulate the central machineries of host cells to their advantage. They prevent host cell antiviral responses to create a favorable environment for their survival and propagation. Measles virus (MV) encodes two nonstructural proteins MV-V and MV-C known to counteract the host interferon response and to regulate cell death pathways. Several molecular mechanisms underlining MV-V regulation of innate immunity and cell death pathways have been proposed, whereas MV-C host-interacting proteins are less studied. We suggest that some cellular factors that are controlled by MV-C protein during viral replication could be components of innate immunity and the cell death pathways. To determine which host factors are targeted by MV-C, we captured both direct and indirect host-interacting proteins of MV-C protein. For this, we used a strategy based on recombinant viruses expressing tagged viral proteins followed by affinity purification and a bottom-up mass spectrometry analysis. From the list of host proteins specifically interacting with MV-C protein in different cell lines, we selected the host targets that belong to immunity and cell death pathways for further validation. Direct protein interaction partners of MV-C were determined by applying protein complementation assay and the bioluminescence resonance energy transfer approach. As a result, we found that MV-C protein specifically interacts with p65-iASPP protein complex that controls both cell death and innate immunity pathways and evaluated the significance of these host factors on virus replication.

Uncovering the genetic blueprint of the C. elegans nervous system.

Despite rapid advances in connectome mapping and neuronal genetics, we lack theoretical and computational tools to unveil, in an experimentally testable fashion, the genetic mechanisms that govern neuronal wiring. Here we introduce a computational framework to link the adjacency matrix of a connectome to the expression patterns of its neurons, helping us uncover a set of genetic rules that govern the interactions between neurons in contact. The method incorporates the biological realities of the system, accounting for noise from data collection limitations, as well as spatial restrictions. The resulting methodology allows us to infer a network of 19 innexin interactions that govern the formation of gap junctions in Caenorhabditis elegans, five of which are already supported by experimental data. As advances in single-cell gene expression profiling increase the accuracy and the coverage of the data, the developed framework will allow researchers to systematically infer experimentally testable connection rules, offering mechanistic predictions for synapse and gap junction formation.

Emergence of disconnected clusters in heterogeneous complex systems.

Percolation theory dictates an intuitive picture depicting correlated regions in complex systems as densely connected clusters. While this picture might be adequate at small scales and apart from criticality, we show that highly correlated sites in complex systems can be inherently disconnected. This finding indicates a counter-intuitive organization of dynamical correlations, where functional similarity decouples from physical connectivity. We illustrate the phenomenon on the example of the disordered contact process (DCP) of infection spreading in heterogeneous systems. We apply numerical simulations and an asymptotically exact renormalization group technique (SDRG) in 1, 2 and 3 dimensional systems as well as in two-dimensional lattices with long-ranged interactions. We conclude that the critical dynamics is well captured by mostly one, highly correlated, but spatially disconnected cluster. Our findings indicate that at criticality the relevant, simultaneously infected sites typically do not directly interact with each other. Due to the similarity of the SDRG equations, our results hold also for the critical behavior of the disordered quantum Ising model, leading to quantum correlated, yet spatially disconnected, magnetic domains.

Scaling of local persistence in the disordered contact process.

We study the time dependence of the local persistence probability during a nonstationary time evolution in the disordered contact process in d=1, 2, and 3 dimensions. We present a method for calculating the persistence with the strong-disorder renormalization group (SDRG) technique, which we then apply at the critical point analytically for d=1 and numerically for d=2,3. According to the results, the average persistence decays at late times as an inverse power of the logarithm of time, with a universal dimension-dependent generalized exponent. For d=1, the distribution of sample-dependent local persistence is shown to be characterized by a universal limit distribution of effective persistence exponents. Using a phenomenological approach of rare-region effects in the active phase, we obtain a nonuniversal algebraic decay of the average persistence for d=1 and enhanced power laws for d>1. As an exception, for randomly diluted lattices, the algebraic decay remains valid for d>1, which is explained by the contribution of dangling ends. Results on the time dependence of average persistence are confirmed by Monte Carlo simulations. We also prove the equivalence of the persistence with a return probability, a valuable tool for the argumentations.

Correlated clusters of closed reaction centers during induction of intact cells of photosynthetic bacteria.

Antenna systems serve to absorb light and to transmit excitation energy to the reaction center (RC) in photosynthetic organisms. As the emitted (bacterio)chlorophyll fluorescence competes with the photochemical utilization of the excitation, the measured fluorescence yield is informed by the migration of the excitation in the antenna. In this work, the fluorescence yield concomitant with the oxidized dimer (P+) of the RC were measured during light excitation (induction) and relaxation (in the dark) for whole cells of photosynthetic bacterium Rhodobacter sphaeroides lacking cytochrome c2 as natural electron donor to P+ (mutant cycA). The relationship between the fluorescence yield and P+ (fraction of closed RC) showed deviations from the standard Joliot-Lavergne-Trissl model: (1) the hyperbola is not symmetric and (2) exhibits hysteresis. These phenomena originate from the difference between the delays of fluorescence relative to P+ kinetics during induction and relaxation, and in structural terms from the non-random distribution of the closed RCs during induction. The experimental findings are supported by Monte Carlo simulations and by results from statistical physics based on random walk approximations of the excitation in the antenna. The applied mathematical treatment demonstrates the generalization of the standard theory and sets the stage for a more adequate description of the long-debated kinetics of fluorescence and of the delicate control and balance between efficient light harvest and photoprotection in photosynthetic organisms.

A reference map of the human binary protein interactome.

Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype-phenotype relationships1,2. Here we present a human 'all-by-all' reference interactome map of human binary protein interactions, or 'HuRI'. With approximately 53,000 protein-protein interactions, HuRI has approximately four times as many such interactions as there are high-quality curated interactions from small-scale studies. The integration of HuRI with genome3, transcriptome4 and proteome5 data enables cellular function to be studied within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying the specific subcellular roles of protein-protein interactions. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms that might underlie tissue-specific phenotypes of Mendelian diseases. HuRI is a systematic proteome-wide reference that links genomic variation to phenotypic outcomes.

The EntOptLayout Cytoscape plug-in for the efficient visualization of major protein complexes in protein-protein interaction and signalling networks.

MOTIVATION: Network visualizations of complex biological datasets usually result in 'hairball' images, which do not discriminate network modules. RESULTS: We present the EntOptLayout Cytoscape plug-in based on a recently developed network representation theory. The plug-in provides an efficient visualization of network modules, which represent major protein complexes in protein-protein interaction and signalling networks. Importantly, the tool gives a quality score of the network visualization by calculating the information loss between the input data and the visual representation showing a 3- to 25-fold improvement over conventional methods. AVAILABILITY AND IMPLEMENTATION: The plug-in (running on Windows, Linux, or Mac OS) and its tutorial (both in written and video forms) can be downloaded freely under the terms of the MIT license from: http://apps.cytoscape.org/apps/entoptlayout. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Publisher Correction: Network-based prediction of drug combinations.

The original version of this Article contained an error in Acknowledgements, which incorrectly omitted the following: 'This work was also supported by NHLBI grant P01HL132825 to A.-L.B.' This has been corrected in both the PDF and HTML versions of the Article.

Network-based prediction of drug combinations.

Drug combinations, offering increased therapeutic efficacy and reduced toxicity, play an important role in treating multiple complex diseases. Yet, our ability to identify and validate effective combinations is limited by a combinatorial explosion, driven by both the large number of drug pairs as well as dosage combinations. Here we propose a network-based methodology to identify clinically efficacious drug combinations for specific diseases. By quantifying the network-based relationship between drug targets and disease proteins in the human protein-protein interactome, we show the existence of six distinct classes of drug-drug-disease combinations. Relying on approved drug combinations for hypertension and cancer, we find that only one of the six classes correlates with therapeutic effects: if the targets of the drugs both hit disease module, but target separate neighborhoods. This finding allows us to identify and validate antihypertensive combinations, offering a generic, powerful network methodology to identify efficacious combination therapies in drug development.

Network-based prediction of protein interactions.

Despite exceptional experimental efforts to map out the human interactome, the continued data incompleteness limits our ability to understand the molecular roots of human disease. Computational tools offer a promising alternative, helping identify biologically significant, yet unmapped protein-protein interactions (PPIs). While link prediction methods connect proteins on the basis of biological or network-based similarity, interacting proteins are not necessarily similar and similar proteins do not necessarily interact. Here, we offer structural and evolutionary evidence that proteins interact not if they are similar to each other, but if one of them is similar to the other's partners. This approach, that mathematically relies on network paths of length three (L3), significantly outperforms all existing link prediction methods. Given its high accuracy, we show that L3 can offer mechanistic insights into disease mechanisms and can complement future experimental efforts to complete the human interactome.